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ABSTRACT 

Understanding atmospheres of protoplanetary cores is crucial for determining the conditions under 
which giant planets can form by nucleated instability. We systematically study quasi-static atmo- 
spheres of accreting protoplanetary cores for different opacity behaviors and realistic planetesimal 
accretion rates in various parts of protoplanetary nebula. We demonstrate that there are two im- 
portant classes of atmospheres: (1) those having outer convective zone which smoothly merges with 
the surrounding nebular gas, and (2) those possessing almost isothermal outer radiative region which 
effectively decouples atmospheric interior from the nebula. A specific type of gaseous envelope which 
accumulates around a given core depends only on the relations between the Bondi radius of the core, 
photon mean free path in the nebular gas, and the luminosity radius (roughly the size of the sphere 
which can radiate accretion luminosity of the core at an effective temperature equal to the local neb- 
ular temperature). Cores in the inner parts of protoplanetary disk (within roughly 0.3 AU from the 
Sun) have large luminosity radii resulting in the atmospheres of the first type, while cores in the 
giant planet region (beyond several AU) have small luminosity radii and always accumulate massive 
atmospheres of the second type. Critical core mass needed for the nucleated instability to commence 
is found to vary considerably as a function of distance from the Sun. This mass is 5 — 20 M^ at 
0.1 — 1 AU which is too large to permit the formation of "hot Jupiters" by nucleated instability near 
the cores that have grown in situ. In the region of giant planets critical core mass depends on the 
gas opacity and planetesimal accretion rate but is insensitive to the nebular temperature or density 
provided that the opacity in the outer radiative region does not depend on the gas density (e.g. dust 
opacity). This is true irrespective of whether the envelope's interior is convective or radiative and 
numerical values of the critical core mass are similar in the two cases. Critical mass in the region of 
giant planets can be as high as 20 — 60 M^ (for opacity 0.1 cm^ g~^) if planetesimal accretion is fast 
enough for protoplanetary cores to form prior to the nebular gas dissipation. This might indicate that 
giant planets in the Solar System have gained massive gaseous atmospheres by nucleated instability 
only after their cores have accumulated most of the mass in solids during the epoch of oligarchic 
growth, subsequent to which planetesimal accretion slowed down and cores became supercritical. 
Subject headings: planets and satellites: formation — solar system: formation — planetary systems: 
protoplanetary disks 



1. INTRODUCTION. 

Recent discoveries of Jupiter-like planets around other 
stars have boosted up efforts to understand the origin 
of giant planets, subject which has been under scrutiny 
since sixties (Safronov 1969; see Brush 1990 for a his- 
torical perspective). Currently one of the most popu- 
lar and successful theories of giant planet formation is a 
core (or nucleated) instability hypothesis (Harris 1978; 
Mizuno et al. 1978). According to this idea massive hy- 
drogen/helium atmospheres of planets like Jupiter and 
Saturn have been acquired as a result of unstable gas ac- 
cretion onto a preexisting core made of rock and/or ice. 
Analytical arguments (Stevenson 1982; Wuchterl 1993) 
and numerical calculations (Mizuno 1980; Ikoma et al. 
2000) suggest that the onset of core instability occurs 
when the mass of the gaseous atmosphere around the 
core becomes comparable to the core mass itself. 

The following gedankenexperiment can help under- 
stand the nature of this instability: imagine placing a 
massive solid core into an infinite homogeneous gaseous 
medium. Gravitational pull from the core gives rise to a 
significant pressure perturbation near the core if the es- 



cape speed from its surface is larger than the gas sound 
speed, see ^2 2.21 On a short (dynamical) timescale gas 
settles into a pressure equilibrium configuration near the 
core. Since the thermal timescale is usually longer than 
the dynamical timescale, this gaseous envelopC"'^ initially 
has entropy equal to that of the gas in the surrounding 
nebula. Gas temperature at the core surface is rather 
high and this drives an outward transport of energy caus- 
ing entropy of the envelope to decrease. As a result, gas 
density near the core goes up and envelope becomes more 
and more massive by slow accumulation (on a thermal 
timescale) of gas from the surrounding nebula. 

Given enough time, atmosphere around the core that 
is not too massive cools down to the nebular temper- 
ature and becomes isothermal, acquiring large gaseous 
mass. Such atmosphere has exponential density profile 
with most of the mass concentrated near the core sur- 
face (Sasaki 1989). Cores are "not too massive" when 
the mass of this isothermal atmosphere is much smaller 
than the core mass. At the same time, exponential sensi- 
tivity to the core gravity (and mass) causes atmospheric 
mass to increase faster than the first power of the core 
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We use terms "envelope" and "atmosphere" intcrchangingly. 
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mass. Because of that, different fate awaits gas around 
the core that is more massive: as the atmosphere cools 
down its mass at some point becomes comparable to the 
core mass and gas starts contributing significantly to the 
planetary gravity. Hydrostatic equilibrium cannot be es- 
tablished beyond this point, because accretion of more 
gas, which helps to reestablish pressure equilibrium, also 
acts now to increase the gravitational acceleration. As 
a result, instability commences allowing gas to accrete 
rapidly onto the core. Critical core mass at which insta- 
bility becomes possible is thus set by the condition that 
the mass of the equilibrium atmosphere around the criti- 
cal core is comparable to the core mass itself [apart from 
additional logarithmic factors intrinsic of the isothermal 
case, see Sasaki (1989)]. Note that if the core is very 
massive, or gas around it is very dense, then even the 
mass of the isentropic envelope forming around the core 
on the dynamical time after the core has been introduced 
into the nebula might exceed the core mass, making at- 
mosphere unstable from the very start. 

Real protoplanetary cores are not likely to be just pas- 
sive solid bodies of constant mass. They accrete planetes- 
imals and accretion not only increases the core mass with 
time but also leads to the release of substantial amount of 
energy at the surface of the core. This changes the steady 
state structure of the envelope which can no longer be 
purely isothermal (as it tends to become in the case of 
passive core). Even if we forget for a moment about 
the increase of the core mass caused by planetesimal ac- 
cretion^, luminosity coming from the core would cause 
temperature in a steady state atmosphere to be higher 
near the core surface than in the nebula. As a result, 
the gas at the same pressure is less dense around the ra- 
diating core than around passive one, meaning that the 
total atmospheric mass is lower for luminous core than 
for passive one. Nevertheless, similar to the isothermal 
case, the more massive the core is, the more massive 
is the envelope around it, and scaling between them is 
faster than linear. As a result, a concept of critical core 
mass holds again: stable equilibrium atmospheres cannot 
exist around the cores which are so massive that their at- 
mospheric mass exceeds the core mass. 

Previous numerical and analytical studies confirm this 
general picture of the core instability, but there are 
still some open issues. For example, Perri & Cameron 
(1974) and later Wuchterl (1993) have concluded that 
protoplanetary atmospheres are likely to be convective 
and that their interior structure and mass (which was 
found to be quite low) sensitively depend on the density 
and temperature of the gas in the surrounding nebula. 
Stevenson (1982) arrived at a very different conclusion by 
considering a simple analytical model of an atmosphere 
with a constant opacity. He found it to be radiative, mas- 
sive, and its structure to be insensitive to the external 
conditions. His findings confirmed previous numerical re- 
sults by Harris (1978), Mizuno ct al. (1978), and Mizuno 
(1980). 

The purpose of this study is to systematically explore 
the structure of atmospheres around accreting protoplan- 
etary cores under a variety of conditions typical in the 

^ Energy release at the core surface might also come from ra- 
dioactive decay or differentiation in the core which are not accom- 
panied by the change in the core mass. 



protoplanetary nebulae. Using analytically tractable but 
still realistic models we attempt to resolve the aforemen- 
tioned issues concerning the state of the envelope. We try 
to single out physical parameters crucial for determining 
the atmospheric structure and this provides us with a 
general classification scheme of possible protoplanetary 
atmospheres. We are looking for the steady state struc- 
tures under the explicit assumption that the mass of the 
gaseous envelope is smaller than the mass of the core 
which it surrounds. This restricts our quantitative re- 
sults to cores with masses below critical mass, although 
all qualitative conclusions are valid even for critical cores 
and we use this to estimate the masses of critical cores, 
see ^5331 

We start by laying down the basics of the problem at 
hand — equations, important length and mass scales, 
boundary conditions — in ^21 In fJ3| & 0| we derive so- 
lutions for two important classes of envelope structures 
typical in the region of giant and terrestrial planets. En- 
velope masses and critical core mass are calculated in Sj^l 
Our results and their implications are discussed in ^ 
Finally, we devote Appendices to technical issues which 
emerge in our calculations. 

2. PROBLEM SETUP. 

Throughout this study the following approximation to 
the protoplanetary disk structure [similar to the Mini- 
mum Mass Solar Nebula (MMSN)] is used: 

Y.g{a) w 100i;p(a) w 3000 g cm^^ „-3/2^ ^^-^ 
To(a) « 300 K a~^/^ (2) 

where Sp,!]^ are the particulate and gas surface den- 
sities correspondingly, Tq is the gas temperature, and 
a,i = a/{n AU) is a distance from the Sun a scaled 
by n AU. From (O one can find the gas sound speed 
Co = {kTo/ jj.y^'^ [k is a Boltzmann constant and is a 
mean molecular weight) and gas density at the midplane 
po = Sgfi/co [SI = (GMo/a^)^/^ is the orbital angular 
frequency, Mq is the Solar mass] : 

Co (a) « 10^ cm s"^ a~^'^, (3) 

pa{a) w 6 X 10-^ g cm^^ (4) 

All numerical estimates in the text refer to this particular 
model of the protoplanetary disk. 

2.1. Basic equations. 

Our purpose is to calculate a spherically symmetric dis- 
tribution of gas density p, pressure P, and temperature 
T around a solid core with the mass embedded in 
the nebular gas. Spherical symmetry requires nebula to 
be at least roughly homogeneous around the core and we 
determine the conditions for this to be true in the next 
section. Structure of the static envelope as a function of 
distance r from the center of the core is governed by the 
equation of hydrostatic equilibrium: 



where G is the gravitational constant and is the mass 
of the core, which we take to be much larger than the 
mass of the gaseous envelope. Equation jSj) is not appli- 
cable if envelope is rapidly rotating (when the azimuthal 
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velocity of gas is of the order of the local Keplerian ve- 
locity around the core), thus in the following discussion 
we assume envelope rotation to be slow. 

Atmosphere is assumed to be heated from below by a 
source at the core surface with luminosity L. This energy 
can be transported by radiative diffusion or convection. 
We use the Schwarzschild criterion to determine the con- 
vective stability of the envelope: 



V < Vad 



7 - 



1 
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ainT 
91nP' 



(6) 



where V is a temperature gradient, Vad is its value under 
isentropic conditions. As usual, 7 is a adiabatic index of 
the gas; for monoatomic gas Vad = 2/5, for diatomic 
(e.g. H2) Varf = 2/7. Usage of the Schwarzschild crite- 
rion implies that envelope is chemically homogeneous and 
nonrotating — a more general Ledoux criterion should 
be used in the presence of the molecular weight gradient 
(Kippenhahn & Weigert 1990), and H0iland criterion has 
to be employed if envelope rotates rapidly (Tassoul 1978). 
When atmosphere is convectively stable according to 
energy released at the core surface has to be car- 
ried away radiatively. In this case one supplements ^ 
with the equation of radiation transfer. In the diffusion 
approximation, valid in the optically thick case, it reads 



IQaT^ dT 
inp dr 



L 



47^7.2 ' 



(7) 



where cr is a Stefan-Boltzmann constant and n is the 
opacity. In the outer Solar System gas around the pro- 
toplanet can be so rarefied that the outer parts of the 
envelope are optically thin. This possibility is consid- 
ered in more detail in H3 3.21 

Whenever the stability criterion © is violated energy 
in the envelope is transported by convection. In this 
study we assume that convection is so effective that the 
temperature gradient in the convective parts of the enve- 
lope is equal to the adiabatic temperature gradient Vad- 
This is equivalent to supplementing (jSJ with adiabatic 
equation of state (isentropic gas) 



P = Kp\ 



(8) 



where K is the adiabatic constant — measure of the gas 
entropy. This approximation should be good enough in 
the dense regions of the envelope (very good in the in- 
teriors of present day giant planets). In Appendix IdI we 
determine under what circumstances this zero entropy 
gradient assumption is valid in the atmospheres of pro- 
toplanetary cores. 

We suppose the luminosity L to be derived from the ac- 
cretion of planetesimals and neglect the additional energy 
release due to the radioactive heating and differentiation 
inside the core. Thus we take 



L = G- 



.MrM 



Rr. 



(9) 



where M is a planetesimal accretion rate. In Appendix 
0we briefly summarize three different regimes of plan- 
etesimal accretion important for the core growth, and 
calculate M and accretion timescale for each of them. 
Equation assumes that (a) planetesimal velocity at 
infinity relative to the core is not too large compared to 
the core's escape speed, and (b) planetesimals penetrate 



to the core surface without much resistance from the en- 
velope and release there all their kinetic energy. First 
assumption is quite reasonable during the buildup of the 
core by planetesimal accretion; second should be valid for 
large planetesimals but small ones may be slowed down 
by the gas drag in the envelope as they make their way 
to the core surface. In the latter case accretion energy 
release does not occur exactly at the core surface but is 
distributed throughout the envelope, meaning that lumi- 
nosity depends on r. However, even in the most unfavor- 
able case of small planetesimals which are quasi-statically 
lowered from the top of the atmosphere to its bottom, 
luminosity is (1 — Rc/r)L, where Rc is the core radius, 
i.e. luminosity is not constant only very near the core's 
surface (Pollack et al. 1996). Beyond several i?c, in the 
bulk of the envelope, we can still safely assume that L is 
constant and given by 

Use of the steady state equations © and lO tac- 
itly assumes that envelope can quickly adjust to the 
changes in the core mass Mc and luminosity L caused by 
the planetesimal accretion. In other words, these equa- 
tions hold only provided that the dynamical and thermal 
timescales of the envelope are shorter than the core ac- 
cretion timescale, and we demonstrate in Appendix [CI 
that this assumption is reasonable. We also show there 
that in the framework of quasi-stationary approximation 
energy release within the envelope due to gas accretion 
can be naturally neglected compared to the planetesimal 
accretion luminosity ©. 

2.2. Important length scales. 

There are several characteristic length scales which are 
important for the problem at hand. One of them is the 
so-called Hill radius r = Rh defined by 

1/3 / n,f \ 1/3 

cm ai 



i?H = a ( ^ I « 2 X 10" cm ai ( ^ ) ,(10) 



Fluid elements at about Rh from the protoplanetary core 
are equally affected by the core gravity and the tidal field 
of the central star. Within the Hill sphere gas dynamics 
is determined by the gravity of the core, while outside of 
it gas is subject mostly to the stellar gravity. 

Physical size of the core Rc scales with Mc in the same 
way as R^ does, meaning that their ratio p is a constant 
depending only on the physical density of the protoplanet 
Pc and the core location in the protoplanetary disk: 



Rc 



/ 3 Mq\ 



-3 



where pi = p/{l g cm^"^). One can easily see that Rc <C 
Rh- 

Bondi radius i?B is defined as the distance from the 
protoplanet at which the thermal energy of the nebu- 
lar gas is of the order of its gravitational energy in the 
potential well of the core: 



Rb — G — ^ 



4 X 10^" cm a\'^ 



Mc 

Mm 



(12) 



Outside Bondi sphere (r > Rb) gravity of the core is 
too weak to strongly affect the gas; consequently, gas 
pressure is almost equal to its nebular value Pq. Inside 
Bondi sphere pressure is significantly perturbed by the 
protoplanetary gravity. 
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Fig. 1. — Different length scales important for atmospheric struc- 
ture as a function of semi- major axis a for kq = 0.1 cm^ i 
and two values of core mass: Ale = 10 {top) and Mc = 1 

Mg) (bottom). Solid, dashed, and dot-dashed lines represent Rg, 
A, and correspondingly. Dotted lines correspond to luminos- 
ity radius Ri^ evaluated for three different planetesimal accretion 
regimes — fast, intermediate, and slow (see Appendix — la- 
beled on the plot. The meaning of other curves is clear from the 
corresponding labels. 



Opacity of the gas sets another important scale. 
We consider a rather general opacity law by assuming 
throughout this work that 



K = ^o {P/Pof {T/ToY 



(13) 



Previous studies (Mizuno 1980; Stevenson 1982) have al- 
ways assumed dust opacity in the outer layers of the enve- 
lope assuming it to be constant. At the same time, opac- 
ity due to small interstellar dust grains (smaller than the 
typical wavelength of the local blackbody radiation) is 
thought to behave as k cx i.e. a = 0) with /3 « 1 — 2 
(Draine 2003) in the range of temperatures typical for 
protoplanetary disks (although this is a statement de- 
pending on the size distribution and composition of the 
dust, which in the protoplanetary disks can be different 
from those in the ISM). Our opacity prescription (|13|l ac- 
counts for this possibility and is more general compared 
to the previous treatments. Opacity in the inner layers of 
the envelope (likely dominated by the molecular opacity 
of H2 and H2O for which one has to use opacity tables) 
might not be so important for the envelope structure or 
mass — atmosphere can be convective there, see t l6 6.11 
Opacity effects are quantified by introducing a mean ab- 
sorption length of photons in the nebular gas 



A = (koPoY 



1.7 X lO-" cm aY 



n/4 / O.lcm^ g-i \ 



Kq 



1(14) 



The exact value of kq is highly uncertain because the 
amount of dust and its size distribution in protoplan- 
etary disks are poorly constrained. Thus, we treat kq 
as a parameter and for simplicity take it to be constant 



throughout the nebula, independent of a. Outer parts 
of the protoplanetary nebula are optically thin, mean- 
ing that A is larger than the vertical disk scale height 
h = co/fl. 

Finally, luminosity of the core sets one more length 
scale 



Rl = 



L 



1/2 



(15) 



which is a size of object radiating luminosity L at an 
effective temperature Tq. In Figuren]we plot i?^ together 
with other length scales as a function of a for two values 
of Kq a-iid for different accretion regimes (see Appendix 

El). 

Nebula can be considered homogeneous on the scale of 
Rb only if Rb is small compared to the disk scaleheight 
h. This condition puts the following constraint on the 
mass of the core: 



12 Me a^/". 



(16) 



Whenever this condition is fulfilled the Bondi radius is 
also smaller than Rh, while Rh is smaller than h, see 
equation (|17|l . For Mc ^ Mi Bondi radius lies outside 
Rh and nebular gas in the Hill sphere is very inhomo- 
geneous {Rh ^ h) and strongl y p erturbed by the pro- 
toplanetary gravity. At 10 AU H16|l constrains Mc to be 
less than « 70Me, while at 30 AU Mi w ISOM®. Note 
that (|16|l is also an approximate condition for the ab- 
sence of strong spatial gradients in the nebula caused by 
the dissipation of the core-induced density waves (Lin & 
Papaloizou 1993; Rafikov 2002). As a result, one does 
not have to worry about the gap formation. 

With the use of H16(l one can rewrite H12() in the fol- 
lowing form: 



Rb 



,'Mc 



2/3 



2/3 



where 



M2 = p^/^Mi « 4.5 X 10" 



-3/4 -1/2 



(17) 



(18) 



is another fiducial mass scale. If the core mass Mc is 
smaller than A/2 then Rb ^ Rc, meaning that the pro- 
toplanet is too small to induce appreciable pressure per- 
turbations in the surrounding gas even at its surface and 
thus has no atmosphere associated with it. We will study 
only the cores satisfying 



M2 < Mc < Ml 



(19) 



since they are massive enough to possess atmospheres 
and small enough for the surrounding gas to be thought 
of as roughly homogeneous on the scale of Rb, i-c. Rc ^ 
Rb <C Rh ^ h. One can see that the mass range in 
which H19I) is valid is rather large and spans 4 — 5 orders 
of magnitude depending on a (e.g. from ~ 0.1 Lunar 
mass to ^ Mj at 10 AU). In Figure El we demonstrate 
how the different length scales are related to each other 
as a function of a and Mc/Mi, constrained by H19(l . 

2.3. Boundary conditions. 

Boundary conditions to the equations of >12 2.11 spec- 
ify that gas pressure, temperature, and density should 
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Fig. 2. — Relationships between different length scales as a func- 
tion of semi-major axis a in the nebula for two values of gas opacity 
(assumed constant throughout the nebula): kq = 0.1 cm^ g~i- (top) 
and kq = 1 cm^ g~^ {bottom). Only the case M2 < Mc < Mi is 
considered (dot-dashed boundaries). Dotted lines correspond to 
condition Rb^ = R\ for three different planetesimal accretion 
regimes labeled on the curves. Solid line represents Rb = A sep- 
aratin g optically thick and optically thin regimes (see H3 3.11 and 
13 3.21 . while dashed line is for Rb = Rl- Squares mark positions 
of 8 major planets of the Solar System on this diagram; for Jupiter 
and Saturn we assume cores of 5 Men and 10 Mm . 



reduce to their nebular values Pq, Tq, and po at some 
distance Rout from the core: 

T{Rout) = To, P{Rout) = Po, p{Rout) = po. (20) 

The problem then is to determine the value of Rout- 

The most simplistic way of evaluating Rout would be to 
completely neglect the fact that the nebula surrounding 
the core is in the differential motion caused by the gravity 
of the central star. If one does this and assumes gas to be 
static everywhere, then Rout can be set equal to infinity. 
Presence of the shearing motion in the nebula drastically 
changes the dynamical and thermal behavior of the gas 
in the core's vicinity. First, gas inside the Hill sphere 
moves in a rather complicated way as demonstrated by 
the the hydrodynamical simulations of planet-disk inter- 
action (e.g. D'Angelo et al. 2002, 2003); but we can hope 
that this will not violate our assumption of nebula homo- 
geneity in the core vicinity if Rb ^ Rh (or Mc ^ Mi). 

Second, the flow of the nebular gas within the Hill 
sphere can effectively "cool" the outer part of the enve- 
lope by advecting the gas heated by the core's accretional 
energy release away from the core and bringing in fresh 
gas having entropy equal to the nebular entropy. This 
process should determine the value of Rout at which tem- 
perature drops to To. Apparently, Rout ^ Rb since gas 
within Bondi sphere is confined by the core's gravity and 
is not being refreshed. At the same time, the distance 
from the core center at which the boundary condition for 
the pressure is imposed is determined by a different pro- 
cess — dominance of the core gravity over the thermal 



pressure in the nebula (and, possibly, over the centrifu- 
gal forces induced by the fluid motion around the grav- 
itationally bound part of the envelope which we neglect 
in this study). Apparently, pressure may converge to Pq 
only beyond Rb', thus, we again arrive at the restriction 
Rout ^ Rb but now from the point of view of pressure 
balance. 

A proper estimate of distance Rout at which the ther- 
mal and pressure boundary conditions are imposed must 
include inhomogeneous heating within Rout, geometry of 
the flow near the core, vertical structure of the nebula, 
possibility of formation of rotationally supported disk 
around the core, and so on. This is beyon d the scope 
of this study. Fortunately, it turns out (see ti6 that 
the exact value of Rout is not important for the envelope 
structure as long as Rout ^ Rb (we checked this to be 
true); in our calculations we simply set Rout ~ 20Rb. 

3. ATMOSPHERES WITH OUTER RADIATIVE ZONE. 

Relationships between the various length scales deflned 
in tl2 2.21 are quite different in the inner and outer parts 
of the nebula. As Figure[5]demonstrates, for protoplane- 
tary cores with Mc « (1 — 1O)M0 accreting planetesimals 
in the fast regime at a > (2 — 5) AU, one of the two sets 
of inequalities 



A ^ Rb, R]^ ^ 
i?L < Pb < A, 



(21) 
(22) 



is typically fulfilled depending on a and kq. This is 
mainly because the planetesimal accretion rate M and 
nebular gas density po decrease with a reducing Rl and 
increasing A. For smaller cores or slower accretion regime 
these conditions can hold even closer to the Su n. We 
demonstrate later on that when either (|22|l or H22|l is 
valid, there exists an almost isothermal radiative layer in 
the outer part of the envelope which decouples its inte- 
rior from the nebular gas. Although the structure of this 
layer is slightly different in two cases, properties of the 
inner envelope will be shown to be universal. 

3.1. Low luminosity, optically thick ('A < i?B and 
Rl « Xtb). 

Whenever the photon mean-free path A is shorter than 
the Bondi radius Rb, envelope is optically thick to the es- 
caping radiation. In this case its structure is completely 
determined by equations lO and lO. Integrating them 
together with (|13|) and using boundary condition P = Pq 
when T = Tq (at r = Rout, see [211]) one finds that 



l+a 



Vo 



- 1 = 



a 



4Vo Rb\ 
3 Rl 



4-f3 



- 1 



4-/3 



(23) 



Constant Vo is positive whenever T increases with depth; 
we assume this to be the case which requires /? < 4 since 
normally opacity does not decrease with increasing den- 
sity (i.e. a > 0). The meaning of Vo will become clear 
later on (see ^'6 3.3|l . The coefficient in the right hand 
side of (|23|) is very large because of H21|) . so that a large 
change of pressure results in only a small perturbation of 
temperature. 
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Fig. 3. — Atmospheric structure around the core with Mc = 
5 at 5 AU. Relative deviations of P, p, and T from their 

nebular values calculated numerically are shown (Pq = 3.2 X 10^^ 
bar, pQ = 7 X 10~^^ g cm~^ and Tq = 130 K in the surrounding 
nebula). Calculation is done for o = 0, /3 = 1, 7 = 7/5 and k = 0.1 
cm^ g~^- In this particular case A < Rg and RgX/R^ « 300, i.o 
envelope i s opti cally thick everywhere and has an outer radiative 
zone, see il3 3.11 Stars mark the position of the outer edge of the 
inner convective region. 



To make further progress we need to substitute H23|l 
into either lO) or (|7I) and find the P and T dependencies 
on r. This can be done numerically for arbitrary a and 
Here we look at the asymptotic behavior of the envelope 
properties in two limiting cases: in the outer atmosphere 
where T — Tq < Tq, and deep in the envelope where 
r>To. 

In the first case envelope is essentially isothermal; solv- 
ing JSJ under this assumption one obtains 



P/Pq = P/Po = exp ( — - 

r Rout 



(24) 



Temperature profile can then be simply obtained from 
(|23) as 



T-Tq 
To 



Rl 



X < exp 



4(1 + 
(1 + a) 



a) Rb^ 

Rb 
r 



Rb 

Ron.t 



- 1 



(25) 



where boundary conditions (|20|) were again used. One 
can see from ()25|) that temperature perturbation at r ~ 
Rb is small because Rb)^/ R\ ^ 1; at the same time gas 
density and pressure within Bondi sphere grow exponen- 
tially, see (|24() . Gas temperature rises to several times 
To only when the distance to the core center becomes as 
small as 

^""^^""HRBXiRiy ^^^^ 

Apparently, Ra ^ Rb although one should keep in mind 



that Ra depends on RbX/R\ only logarithmically. Pres- 
sure Pa and density pa at this depth are 

PalPo ~ Pa/Po ~ (i?BA/i?i)l/(l+") » 1. (27) 

Existence of the outer radiative layer, in which temper- 
ature is almost constant while density dramatically in- 
creases was first found numerically by Harris (1978) un- 
der the conditions typical in the region of giant plan- 
ets. Note that the presence of the outer isothermal layer 
does not necessarily require gas to be optically thin, cf. 
Mizuno et al. (1978). 

In a second limiting case, interior to Ra, we use T 3> To 
to integrate with Using T « ^Tq (^ 1) at 

r — Ra as an approximate boundary condition at the 
transition point one finds that 



To ^ 



Vn 



P 

Po 



j^ 4Vo XRb 



Rb Rb 

1/(1+q) 



Rb 
r 



Rb 

Ra 



(28) 

l/Vo 

(29) 



In the case of constant opacity (a = /? = 0, Vq = 1/4) 
this solution reduces to the radiative zero solution found 
by Stevenson (1982) for r <C Ra- Figure El shows the 
internal structure of the atmosphere around 5 core 
at 5 AU (for which [2] is valid) calculated using and 
H23() . Opacity with a = 0, (3=1 and kq = 0.1 cm^ g~^ as 
well as 7 = 7/5 are assumed in calculation. This specific 
scaling of k with P and T makes the inner part of the 
envelope convective below Ra (see H3 3.3() , but this hardly 
changes the overall picture of the atmospheric structure 
that we described in this section. 

3.2. Low luminosity, optically thin (Rl <C Rb ^ Aj. 

Whenever A > Rb there is an optically thin region 
around the protoplanet in which the photon mean free 
path is larger than the length scale over which physical 
variables such as pressure and temperature experience 
significant changes. Deep in the envelope (below the 
photosphere) optical depth increases and atmosphere be- 
comes optically thick. Gas is optically thick also far from 
the core since the density and temperature scale height 
is ~ r there and the optical depth to radiation escaping 
from the core becomes larger than unity at r > A. Of 
course, this outer optically thick region can only exist if 
A < ft. and nebula is homogeneous on scales A, which 
is not the case in the outer parts of protoplanetary disks 
(but this turns out not to be important). Thus, the op- 
tically thin zone around the core in the infinite medium 
should be sandwiched between the inner and outer opti- 
cally thick regions. 

Temperature structure in the optically thick parts is 
determined by equation In the optically thin region 
we have T"* w Tq -f L/(167r(Tr^), where the additional 
factor of 4 in the denominator of the second term comes 
from the anisotropy of the core radiation. Similar be- 
havior of T in the optically thin region was found by 
Hayashi et al. (1979). Strictly speaking, this expression 
is accurate only far from the photosphere and we should 
expect it to reproduce photosphcric temperature only ap- 
proximately (which will not affect our results in any sig- 
nificant way). Thus, we express envelope temperature 
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Fig. 4.— Same as Figure El but for Mc = 10 M® at 30 AU. In 
this particular case Pq = 10~^ bar, po = 5 X 10"^'^ g cm~'^ and 
To = 55 K. Also A > Rb and RBX/Rj^ ~ 10^, i.e envelope has an 
out er opti cally thin region and possesses an outer radiative zone, 
see i |3 3.21 The rest of notation is as in Figure 



profile above the photosphere as 

L 3L 



Kp{r')dr' 



(30) 



The second term on the right hand side is relevant for the 
optically thin region while the third becomes important 
outside of r sa A where photons couple to the nebular gas 
and start to diffuse. 

Gas pressure is close to Pq for r > i?^. Assuming that 
T — To <C To everywhere in this region we have p ~ po 
and K w Ko- Then it follows from (|30|) that 

1/4 



T(r) » To 



l + -^+3— ^ 



(31) 



It is very important that the temperature perturbation at 
the Bondi sphere T{Rb)-Tq ~ ToiRL/Rs)^ is negligible 
whenever H22|l is fulfilled (similar to the optically thick 
case considered in ^3 3.1|l — this verifies our assumption 
of r « To for r > Rb- This result is independent of 
whether A is larger or smaller than h or Rjj, what is cru- 
cial is that Rl <C Rb ^ ^ (which makes third term in 
|ST] negligible compared to the second at r ~ i?^). Tem- 
perature deviation T — Tq is dominated by the radiative 
diffusion for r > X (third term in brackets) and by the 
optically thin radiation transfer for r < A. 

Inside the Bondi sphere gas is initially still optically 
thin and its temperature is determined by (|31|l meaning 
that T m Tq. At the same time, pressure and density in 
this essentially isothermal region increase exponentially 
with depth in accordance with 124|) . As a result, local 
photon mean free path rapidly decreases and envelope 
finally becomes optically thick to the escaping radiation. 
Photosphere is located at the distance Rph from the core 



center where the photon mean free path becomes com- 
parable to the typical length scale of density variation. 
Using (|24|) we estimate that the density scale height is 
dr/dlnp — /Rb, which becomes comparable to (kp)~^ 
at 

^'"'"'^^InW^s)' ^^^^ 

whenever Rb <C A. We find from (|31|l that gas temper- 
ature at the photosphere Tph = T[Rph) is offset from Tq 

by 

Tph - To _ Rl _ [In (X/Rb) Rl ' ' 



To 



2(1 -ha) Rb 



< 1, (33) 



while the photospheric pressure Ppu = P{Rph) is 

PphlPo ~ (A/i?B)i/(i+") » 1. (34) 

We see that Tph is still only slightly different from To, 
while the photospheric pressure is much higher than Pq 
if ^ holds. 

Below the photosphere, we again resort to equations 
Q and 10. Similar to (|23|l one finds for boundary con- 
ditions P = Pph dXT — Tph that 



p \ i+a 

4Vo RbX 
3 Rl 



Pph 
Po 
T 



Tph 



4-/3- 



(35) 



When I T- To I /To < 1 we can still use ^ for the density 
profile; then the temperature profile for r < Rph derived 
from H35|) becomes 



T-T, 



ph 



Rl 



exp 



4(1 -ha) RbX 
Rb 



(1 + a) 



{PvH/P,r 



(36) 



where we took into account that Rb/ Rout ^5 1- Keeping 
in mind that according to (|22Jl and (jSUl {Pph/PoY^" < 
RbX/R\, one easily finds that the distance Ra at which 
gas temperature starts to appreciably deviate from To is 
still given by (|26|) . Thus, under conditions 1)22(1 Ra is still 
located deep inside the envelope, and always below the 
photosphere, compare with (|32|l . 

It is also easy to deduce from (|35|) that for r < Ra tem- 
perature and pressure behavior in the inner envelope are 
still described by equations and (j^ . Thus, despite 
the differences in the structure of the outer atmosphere 
{r > Ra) in the optically thin (i?^ <C i?_B ^ A) and opti- 
cally thick (A <C Xtb) cases, the structure of 
the inner atmosphere is the same. In Figure ^properties 
of the envelope around 10 core at 30 AU (for which 
|22) is valid) are exhibited. Temperature profile in the 
optically thin region was calculated using ((30|l . and all 
parameters (opacity, value of 7) are the same as those 
used for Figure |21 see i i3 3.11 

3.3. Convective stability. 

In M3 3.113 tOl we have assumed that energy transfer 
from the bottom to the top of the atmosphere occurs by 
radiative transport. Here we explicitly determine under 
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which circumstances this is the case, and also calculate 
envelope structure in convectively unstable regions. 

To check our solutions obtained in previous sections 
for convective stability we calculate V and use the 
Schwarzschild criterion It turns out that the outer 
part of the envelope above Ra is stable whenever (j21|l or 
(|22|) is fulfilled. Indeed, in the optically thick case (see 
H3 3.1|l we find using I©, 10, and (EH) that 



v(r) 



3 Rl 



ARbX 

To 
T 



T 

4-/3 



13-4. 



1 - 



1+a 



4Vo Rb\ 



(37) 



Because of the limitations imposed by H21f) this reduces 
to 

1 



Vo 



(38) 



[relative corrections to this expression are at the most 
0{RI/RbX) < 1]. Thus, outside of Ra, where T « Tq, 
temperature gradient is small and atmosphere is con- 
vectively stable. We also see from (|38|) that Vq has a 
meaning of temperature gradient deep inside the radia- 
tive envelope where T ^ Tq. 

In the optically thin case (see H3 3.2() one finds similar 
situation. Temperature above the photosphere is almost 
equal to Tq (see equations and while pressure 
profile is given by (|24|) . Then one obtains that 



4 V Rsr 



3«i 



RbXJ 



« 1, 



(39) 



i.e. atmosphere is convectively stable above the photo- 
sphere. Below the photosphere, for r < Rp^, tempera- 
ture is related to pressure by (|35|l which results in 



v(r) = Vo <^ 1 - 



1 



Rl 



4Vo RbX 



Tph 

Pph 

Po 



4-f3 



Tph 



4-/3- 



(40) 



Since according to k Tph « Tq and 

(Pph/PnY^"' ~ X/Rb, one finds that V below Rph is still 
given by H38|l with the relative accuracy 0{R\/ R\) <^ 1. 
Consequently, in the optically thin case (|22|) atmosphere 
is convectively stable everywhere above Ra, analogous to 
the optically thick case. Thus, in both cases considered 
in i j3 3.11 and 13 3.21 there is an outer radiative ( convec- 
tively stable) region in the atmosphere, which is almost 
isothermal. 

Deep in the envelope, below Ra, radiative tempera- 
ture gradient is given by l|38() and steadily increases with 
depth (since T monotonically goes up); deep inside the 
envelope V becomes equal to Vq- Thus, whenever 



Vo 



1 



a ^ 7—1 

< Vad 



1 



(41) 



the whole envelope is convectively stable and energy is 
transferred from the bottom by radiative diffusion. If 
opacity is independent of P and T, i.e. a — (3 — Q, then 
Vo = 1/4 (Stevenson 1982) and envelope is convectively 



stable provided that it is composed of monoatomic or 
diatomic gas. 

In the opposite case, when Vq > Vad envelope becomes 
convectively unstable at some point. From (|38|l we find 
that this happens when gas temperature reaches 



Tn 1 



Vad 
Vo 



-1/(4-/3) 



(42) 



This critical temperature above which convection sets 
in is clearly not very different from Tq. Thus, Tconv is 
achieved at « Ra from the core center and we conclude 
that atmosphere becomes convective below « Ra- For 
example, when the dominant source of opacity in the 
outer atmosphere is dust with f3 — 1, one should ex- 
pect envelope to be convectively unstable if II2 is its 
major constituent. In this case convection sets in at 
Tconv = 71/37^0 ~ i.9To. The edge of the convection zone 
as shown in Figures 13 and 01 agrees with this estimate. 

If envelope is convectively unstable below Ra , equation 
cannot be used there. Instead, one has to utilize 
to relate pressure and density in 0. Value of if in JSJ) is 
set by the conditions at the edge of the convection zone, 
i.e at r = Tconv In the optically thin atmospheres (see 
i i3 3.2|l pressure at this point is set by (P5|l to be 




4Vo RbX 



Rl 



1 - Vad/Vo 



Tph 
To 



4-/3- 



1/(1+") 



Po 



4 VadVo RbX 

3 Vo — Vad Rl 



1/(1+") 



(43) 



The last line in follows from and Tph ~ To. Ap- 
parently, Pconv Pa ^ Po- In the optically thick atmo- 
spheres pressure is given by (|23() , which can be obtained 
from IpT^ by setting Pph = Po and Tph = To. Making the 
same substitutions in equation (|43|l we find that Pconv is 
still given by H43() , despite the different structure of the 
outer radiative layer. 

Solving equations and (0 with the boundary con- 
dition p = Pconv = Pconv lJ-/kTconv at T « i?a and using 



K^P 



p{r) 



'conv Pconv OUC fiuds 



1 + Vad' 



GMc 



7-1 

conv 



1 + V„d 1 - 



For r < Rail - i?a/i?s) 
Pconv and one obtains 



1 1 

r^R^ 

1/(4-/3) 



1/(7-1) 



Rb 
r 



Rb 

Ra 



H 1/(7-1) 

(44) 



gas density strongly exceeds 



Pc 



7-1 



T 



Tc, 



"^ad 



To 
T 



Rb 
r 



Rb 

Ra 



(45) 



Note that this temperature profile would be identical to 
H28I) found for envelopes with radiative interiors if V^d 
were replaced with Vq; apparently, temperature behavior 
is not very sensitive to the details of physics determining 
the atmospheric structure. 
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4. ATMOSPHERES WITH OUTER CONVECTIVE ZONE. 

Whenever planetesimal accretion luminosity is high, 
luminosity radius Rl can exceed Rb in the optically thin 
case ( c.f. ^TO) or {RBXf^ in the optically thick (c.f. 
ii3 Based on the results of SjSlwe can anticipate that 
intense energy release at the core surface would strongly 
affect gas temperature even beyond the Bondi radius, 
and this qualitatively changes the atmospheric structure. 

We will concentrate on the optically thick case with 



Re S> a, i?!^ > Rb^ 



(46) 



(directly opposite to |23), typical in the terrestrial re- 
gion (see Figure 121). Since atmosphere is optically thick 
equation (|23|l should determine the envelope structure if 
it were convectively stable. However, it is easy to show 
that under the conditions H46|) gas around the core can- 
not be convectively stable even at Rout ^ Rb^- Indeed, 
temperature gradient is given by 1)37(1 and one finds that 
far from the core V > Varf provided that 



Vad < V(oo) 



3_R^ 

ARbX 



3 LkqPqcI 
647r GMcCrTfj ' 



(47) 



If gnj holds, V « R\/{Rb\) > 1 already at r - Rout, 
meaning that even the outer part of the envelope is con- 
vective. This is completely different from the situation 
typical for the region of giant planets (see ^ where con- 
ditions are such that protoplanetary envelopes always 
have an almost isothermal outer radiative zone. The 
reason for this difference is that in the present case large 
energy flux escaping from the core severely affects gas 
temperature outside Bondi sphere where gas pressure is 
still almost equal to its nebular value Pq. Consequently, 
temperature gradient d\nT / d\nP takes on a high value 
outside of Rb (unlike the case studied in Q giving rise 
to convection. As a result, envelope acquires an outer 
convective zone. 

Calculation similar to that leading to (|44|l yields the 
following density structure in the convective part of at- 
mosphere: 



P{r) = Po 



Rb Rb 



Rou 



1 1/(7-1) 



(48) 



This expression shows that gas density, pressure, and 
temperature deviate weakly from their nebular values po, 
Pq , and To as long as Rb ^ i" < Rout , but they increase as 
power laws of Rb /r inside the Bondi sphere. Apparently, 
internal atmospheric structure (for r < Rb) is rather 
insensitive to a particular choice of radius Rout at which 
boundary conditions are set, as long as Rout ^ Rb- 
Structure of the interior regions of the envelope de- 
pends on the opacity behavior as described in H3 3.31 at- 
mosphere is convective from Rout all the way to Rc when- 
ever Vo > Vad, but it becomes convectively stable and 
radiative at some depth if Vq < Va^. In the latter case, 
for the inner radiative region to exist it is necessary that 
convection stops above the core surface. Suppose that 
Vo < Vad and that transition between the outer convec- 
tive zone and the inner radiative region takes place at a 
distance Rrad from the core center. Pressure, tempera- 
ture, and density at this point are Prad, Trad, Prad', they 
can be uniquely related to Rrad by H48|) since the outer 
boundary of the radiative zone is also an inner boundary 
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Fig. 5. — Same as Figure|2lbut for M, 

particular case Pq = 3 X 10~^ bar, po = 

To = 550 K. Also X ^ Rb and Rb\/R\ ^ 0.06, i.e envelope is 
optically thick everywhere and possesses an outer convective zone, 
see Jl] The interior of this atmosphere is also convective for a 
particular choice of parameters (the same as for Figure |3 used in 
producing this Figure. 



of the convective zone. Within the radiative zone behav- 
ior of V(T) can be described by H40|l with Pph and Tph 
replaced by Prad and Trad- Then one can easily fix the 
value of Rrad from the condition V(Trad) — ^ ad'- 



Rrad ~ Rb 



e ' 



6; 



/ 3 



Rl 



(<l-|3)(V„^-Vo) 



V4Vad RbXJ 



(.49) 



It is clear that 6 1 and Rrad ^ Rb because of 
14(j|) . Inner radiative zone exists only if Rrad > Rc 
in addition to Vo < Vad- Temperature and density 
at Rrad are, of course, much larger than To and Po: 

Trad/To = (prad/poV-^ ~ &, See 

We also briefly discuss the atmospheric structure in the 
optically thin case'^ 



Rb <C a, Rb <C Rl 



(50) 



opposite to that considered in ^3 3.21 Using equation H31|l 
to determine the temperature structure in the optically 
thin part of the envelope one finds T to be strongly per- 
turbed beyond the Bondi radius because of the second 
inequality in (|50|) . Then it is rather clear that such at- 
mospheres should possess outer convective zones, similar 
to the case when (|46ll holds. This conclusion is indepen- 
dent upon the relationship between R^ and A. We do 
not consider the interior structure of such envelopes in 
this study. 

^ This particular relationship between Rg, ij^, and A can be 
realized only in the region of giant planets for small cores accreting 
in the fast regime, see Figure |2] We mention it here mainly for 
completeness. 
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5. ENVELOPE MASS AND CRITICAL CORE MASS. 

We define mass of tlie envelope M^nv as 



Menv = ^7r j p{r')r"dr', 

Rc 



(51) 



wfiere we have chosen Rb to be the outer boundary of the 
envelope. We compute Me„„ separately for atmospheres 
with the outer radiative and convective zones since it will 
turn out that masses are very different in the two cases. 
Using these results we also estimate the critical core mass 
necessary for the initiation of a runaway gas accretion in 

5.1. Envelopes having outer radiative zone. 

Results of inidemonstrate that gas density in the atmo- 
spheres having outer radiative zone increases exponen- 
tially between Rb and Ra- As a result, most of the at- 
mospheric mass is contained within Ra- In Appendix 
El we demonstrate that Menv is given by 



47r*ipoi?| 



/ Rb^\ 



1/(1+") 



(52) 



where is a weak (logarithmic) function of Rb)^/R\ 
given by (|ljl() . Envelope mass is dominated by the con- 
tribution coming from r ~ R^ in all dynamically sta- 
ble atmospheres with convective interiors and in all at- 
mospheres with radiative interiors having Vq > 1/4; in 
these cases the inner part of the atmosphere near the core 
contributes to Menv only weakly, see Appendix IbI We 
restrict ourselves to studying only these two important 
classes of envelopes. 

Of special interest is the case when atmospheric opacity 
is independent of pressure (and, consequently, density), 
i.e. a = 0. Then, regardless of whether envelope interior 
is radiativ e or convective, one finds using ((T^ . ifTHl . (fT5|l . 
(f^ . and lEU that 
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Different numerical estimates are for three accretion 
regimes described in Appendix 1X1 and assume molecular 
gas of cosmic composition, /? = 1, and 7 = 7/5 (con- 
vective interior). The peculiarity of this special case is 
that envelope mass is virtually independent of the tem- 
perature and density in the surrounding nebula. Indeed, 
both To and po enter H53|) only logarithmically through 
the dependence of on Rs/Ra, sec and (|Blll . Lo- 
cal conditions in the protoplanetary disk affect Menv only 
through gas opacity kq and planetesimal accretion lumi- 
nosity L. This was first noticed by Stevenson (1982) who 
discovered this feature while studying radiative envelopes 
with constant gas opacity, a — f3 — 0. But equation 
demonstrates that insensitivity of Menv to Tq and po is 
a more general phenomenon since it holds even when at- 
mospheric opacity varies with temperature and for con- 
vective as well as radiative interiors, provided only that 
a = 0, i.e. opacity is independent of the gas density. 




Fig. 6. — Mass of the envelope for two values of Mcore = 
1, 10 Afgj as a function of semi-major axis a in the nebula. Different 
panels correspond to different planetesimal accretion regimes: (a) 
slow, (b) intermediate, and (c) fast. Solid curves labeled with cor- 
responding core masses represent results of numerical calculation 
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same given by 1531 and I55i for large and small a respectively. 



When K does depend on pressure {a ^ 0), one finds 
that 



Me. 



PZiMcp) 



4+3a 



1/(1+") 



(54) 



i.e. envelope mass depends on nebular properties in this 
more general case. But even then the detailed character 
of dependence is not determined by whether atmospheric 
interior is radiative or convective, but only by the opacity 
dependence on P. 

It might seem surprising that mass of the envelope with 
convective interior can be sensitive to the gas opacity 
since the energy transfer below Ra is not done by radi- 
ation. The explanation lies in the presence of the outer 
radiative zone above Ra which is non-adiabatic. Because 
of that gas entropy at the edge of the convective zone 
(at « Ra) is set by the radiative energy transfer in the 
outer atmosphere and depends on the gas opacity kq (see 
equation [53 )• This is the origin of dependence of Menv 
on the nebular opacity in the case of atmospheres with 
convective interiors. 

In Figure we plot Menv for cores of specified mass (1 
and 10 Mq) as a function of core's distance from the Sun 
a for fast, intermediate, and slow planetesimal accretion 
regimes (see Appendix IXl . Numerical results of the en- 
velope structure calculations for a = 0,/3 = 1, kq = 0.1 
cm^ g^i, 7 ~ 7/5 are shown by solid curves; dashed 
line at large a is our analytical estimate of Menv given 
by (jnSl- As Figure demonstrates, formula (|53|l some- 
what overestimates Menv'- because of t he fi nite size of 
the core our extension of integration in |IB1|) to zero (in- 
stead of Re) leads to the overestimate of Menv', relative 
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correction is '-^ {Rd RaY^"^ which is smah but sometimes 
non-negUgible. We checked that this discrepancy goes 
away when we artificiaUy set Rc to in our numerical 
calculations. 

5.2. Envelopes having outer convective zone. 

Mass of the atmosphere possessing outer convective 
zone is calculated using density profile (|48|) . Assuming 
that (|48|) holds up to the core's surface (fully convective 
envelope, Vq > Vad) one finds 

3 



Menv = 4^^'2P0i?B « 4 X 10"* g — ^ 



Ma 



-5/4 



*2(7) 



Vad(7) 



1/(7-1) 



(55) 



where we have assumed Rs/Rout S 1- Numerical esti- 
mate is done for 7 = 7/5 (^'2 ~ 0.88). Similar estimate 
of Menv can be found in Wuchtcrl (1993). 

A remarkable property of atmospheres having outer 
convective zone is that gaseous mass contained within 
a Bondi sphere around the core is of the order of the 
mass of nebular gas with p — po contained within the 
same volume! This is very different from atmospheres 
possessing outer radiative region, which not only have 
envelope mass much higher than {A/3)TTpaR'g, but also 
contain this mass within smaller volume than that of the 
Bondi sphere, see (|52|) . In addition, M^nv given by H55|l 
is completely independent of the core luminosity and gas 
opacity, very much unlike the case studied in ijb 5.11 

Most of the mass in fully convective atmospheres is 
concentrated in their outer part (because of the require- 
ment 7 > 4/3 necessary for dynamical stability), which 
allowed us to set the lower integration limit in the defini- 
tion of to zero. One can demonstrate that in the case 
of envelopes having inner radiative region (those with 
Vo > Vad) formula continues to correctly describe 
the mass of the envelope, provided that Vq > 1/4 in its 
radiative part (because most of the mass in radiative re- 
gion is then concentrated near its outer edge); also, inner 
radiative region typically has rather small radial extent, 
see 

In Figureinianalytical estimate (|55|l is displayed by the 
dashed line in the inner part of the nebula (terrestrial 
planet region, < 1 AU). One can see very good agreement 
of analytical result with the predictions of more detailed 
numerical calculations. 

5.3. Critical core mass. 

As we have mentioned in the Introduction, both nu- 
merical calculations and analytical arguments suggest 
that phase of rapid gas accretion onto the protoplane- 
tary core initiates when M^nv ~ Mc. The exact ratio of 
the two masses at the onset of instability is uncertain and 
can be determined only after the envelope's self-gravity is 
self-consistently taken into account which is beyond the 
scope of this study. Here we simply assume this critical 
ratio to be a free parameter ~ 1, so that 

Menv{Mcr) = VMcr, (56) 

where Mcr is the critical core mass at the onset of rapid 
gas accretion. For a fixed rj instability condition (|56|) can 
be viewed as an equation for Mcr. 




Fig. 7. — Critical core mass as a function of semi-major axis a in 
the nebula for (a) slow, (b) intermediate, and (c) fast planctesimal 
accretion rates. Solid curve represents numeric al r esults , w hile 
dashed lines display estimates of Mcr given by 1571 and 1581 at 
large and small a correspondi ngly. Dot-dashed line shows the run 
of fiducial mass Mi (equation 1161 ) with a. 



Using (O, (|53|l . and (|A1|I we find that for atmospheres 
having outer radiative region (Q critical core mass is 
given by 
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where 770.3 = 77/O.3 and is a parameter determining the 
efficiency of accretion. Apparently, core instability sets 
in at lower M,. in the more distant parts of the nebula 
because of the rapid decrease of planctesimal accretion 
rate M with a. Note a strong dependence of Mcr on the 
mean molecular weight p, (Stevenson 1982). Convective 
erosion of the core and dissolution of infalling planetesi- 
mals in the envelope might increase p and considerably 
lower Mcr^ facilitating rapid gas accretion. 

Opacity determining Mcnv and Merit in l|52(l and H57|l 
is the opacity in the outer radiative zone only — it is the 
radiation transfer in this region that determines the con- 
ditions in the innermost part of the envelope and, hence, 
its mass. Temperature in the radiative zone does not de- 
viate strongly from Tq and we can be more confident that 
the opacity behavior in this region can be represented by 
a simple power law dependence ifT^ . This assertion be- 
comes even more robust if opacity there is also indepen- 
dent of gas density (since density strongly varies within 
the radiative zone), and this is the case when opacity is 
dominated by dust absorption. 

In the case of atmosphere having outer convective zone 
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(2J we find after substituting (|55|l into (|56|l that 
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This critical mass decreases as one goes inward in the 
nebula (unlike [S7|) because gas density in the disk 
rapidly increases inward. Introducing the Toomre stabil- 
ity parameter Q = flcg / {ttGTj g) one can rewrite H58|l as 
Mcr = (Q/4^'2)^/^Mi, see definition This impHes 

that Mcr ^ Ml in gravitationally stable protoplanetary 
disks^ with Q > 1. It is then clear that M^nv given by 
(|58l) is only a rough estimate of the critical core mass 
because nebula cannot be considered static and homoge- 
neous on the scale of i?_B (as we always assumed) when 
Mc > Ml, see ^TT^ 

Comparing H52() with (|55|l we see that presence of the 
outer radiative zone, requiring Rb^ 3> R^, increases 
Menv and lowers Mcr- This happens because radiative 
diffusion across almost isothermal outer radiative zone 
significantly reduces the entropy of the inner envelope. 
Lower entropy means higher density resulting in a higher 
Menv compared to the mass which an isentropic atmo- 
sphere would have for the same core mass. For instance, 
in the optically thick atmosphere with RbX ^ R\ the 
value of adiabatic constant K at Ra is 



K 



(7-l)/(l+a) 



(59) 



see 122], 123, and [221; i^o = -PoPo^ is the adiabatic 
constant of gas in the surrounding nebula. It follows 
from (|59|) that K ^ Kq, i.e. gas entropy in the inner 
atmosphere is much lower than it is in the surrounding 
nebula. Consequently, Mcr for an atmosphere with the 
outer radiative zone is lower than it would have been 
if atmosphere were isentropic (i.e. had outer convective 
zone). 

In Figure we plot the critical core mass as a func- 
tion of a for different planetesimal accretion regimes, self- 
consistently taking into account transition between dif- 
ferent types of atmospheres in the inner and outer parts 
of protoplanetary disk. In addition to the results of nu- 
merical calculations of the atmospheric structure we also 
plot analytical approximations for Mcr given by (|58|l in 
the region of terrestrial planets and by H57(l in the region 
of giant planets. In the terrestrial planet region theory 
and numerical calculations agree very well and one can 
see that Mc > Mi in agreement with what we have said 
before (Q > 30 at a < 1 AU for the MMSN parame- 
ters given by 121). In the giant planet region there are 
discrepancies between the theory and numerical results 
at large a which are especially pronounced in Figure C^. 
They appear because of the finite size of the core which 
causes (j5^ to overestimate Mgnv , see MS 5.11 in the more 
careful numerical calculation a bigger core is required to 
trigger the runaway gas accretion. 

6. DISCUSSION. 

One of the major results of our study is a clear distinc- 
tion between protoplanetary atmospheres having convec- 
tive outer region which smoothly merges with the sur- 
rounding nebula and atmospheres having radiative zone 

* Similar argument was advanced in Ikoma et al. (2001). 



separating dense and hot interior from the nebula. A spe- 
cific type of atmosphere around a particular protoplan- 
etary core is determined by the relationships between 
three important length scales — Rl, and A — pro- 
vided that core mass satisfies the mass constraint (I19|l . 
Whenever planetesimal accretion rate of the core is low, 
meaning that Rl is small (i.e. either [22] or [22] is ful- 
filled), gas temperature can be appreciably affected only 
deep inside the Bondi sphere while the pressure starts to 
vary already at the Bondi radius. As a result, according 
to the Schwarzschild criterion ||BJ), outer parts of the en- 
velope are convectively stable and energy is carried away 
by radiation. In the opposite case of very high accretion 
rate and large Rl (i-e. either 0S] or [SO] holds), gas tem- 
perature is perturbed even outside the Bondi sphere, and 
condition © predicts that gas is convectively unstable 
for r > Rb, where pressure perturbation is small. More- 
over, in the latter case even if the innermost parts of the 
envelope tend to settle onto the radiative (convectively 
stable) configuration, they would be able to switch to 
a radiative solution only very deep in the envelope, see 
equation (|^ . 

This segregation of atmospheres into two major classes 
depends to some extent on whether the outer parts of the 
atmosphere are optically thick or thin. For example, en- 
velope structure in the optically thick case depends only 
on the value of the dimensionless parameter XRb/R\ — 
roughly the inverse of the radiative temperature gradi- 
ent far from the core (see ^2], 02], CH], and |D1| '). and 
is completely insensitive to the individual relationships 
between A, Rl, and Rb (as long as A <C Rb)- enve- 
lope has outer convective zone when \Rb/R\ < 1 and 
outer radiative zone when \Rb/R\ > 1. Analogous con- 
dition was formulated by Wuchterl (1993) in terms of 
nebular density po- Separation between the two classes 
of atmospheres in the optically thin case is somewhat 
more complicated (e.g. see constraint [22]) and one has 
to pay attention to the individual relationships between 
the three length scales. Results of |J2]&E] cover all such 
possibilities. 

Because these characteristic length scales vary with 
the distance from the Sun, a specific type of atmosphere 
forming around the core of a given mass depends on a. 
This is caused primarily by the variation of the plan- 
etesimal accretion rate throughout the protoplanetary 
disk: in the terrestrial region the planetesimal surface 
density is high while the dynamical timescale is short 
leading to high M, large Rl, and small Rb^/R?l- An- 
other important factor is a steep dependence of A on a 
(see Figure nj. As a result, atmospheres around massive 
cores in the inner parts of protoplanetary disks (within 
roughly 0.5 — 2 AU depending on the planetesimal accre- 
tion regime) where i^sA < R\ possess outer convective 
zones, see Figure[2t- On the contrary, in the region of gi- 
ant planets both Sp and SI are small, which decreases M 
making Rl small as well. Photon mean free path there 
is large because gas density is very low (disk can be op- 
tically thin). These factors conspire to rapidly increase 
Rb^I R\ with a and allow cores outside of ~ (0.5 — 2) AU 
to have massive atmospheres with quite extended outer 
radiative zones, see Figures [S] & E] Some exceptions are 
possible (e.g. cores having Rb "C A and Rb ^ Rl in 
the region of giant planets have convective outer atmo- 
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spheres, see 21 ^-nd Figure |2) but they typicaUy occur 
for rather small cores unable of retaining massive atmo- 
spheres. 

6.1. Comparison with previous studies. 

Envelopes with the outer convective region have been 
previously considered by Perri & Cameron (1974), 
Wuchterl (1993), and Ikoma et al. (2001), who found the 
masses of such envelopes to be rather small which trans- 
lated into large critical core mass. In addition, prop- 
erties of such fully convective atmospheres were found 
to strongly depend on the temperature and density in 
the surrounding nebula. This is exactly the picture that 
we described in 21 Wuchterl (1993) suggested that the 
sensitivity to external conditions is a consequence of en- 
velopes being fully convective, but as we demonstrate in 
f^lwhat is really important is the convection in the outer- 
most region of atmosphere only, independent of whether 
the atmospheric interior is convective or radiative. Pres- 
ence of the outer convective region sets entropy of the 
inner atmosphere equal to the entropy of the nebular 
gas, so that variations in the external conditions directly 
affect the overall structure (and mass) of the envelope. 
Nebular entropy is quite high, which makes atmospheres 
not very dense and accounts for the low masses of such 
envelopes, see Figure El As our analysis demonstrates, 
Perri & Cameron (1974) and Wuchterl (1993) have prob- 
ably stretched their assumptions too far by calculating 
Mcr ~ 60 for convective envelopes at 5 AU, in the 
region of giant planets, where atmospheres should in fact 
have outer radiative (not convective) zones. 

Envelopes possessing a radiative region between the 
inner, dense parts of the atmosphere and the nebular 
gas outside have first been studied numerically by Harris 
(1978), Mizuno et al. (1978), and Hayashi et al. (1979). 
Mizuno (1980) and Stevenson (1982) were the first to 
notice that presence of the outer radiative region makes 
critical core mass virtually independent upon the den- 
sity and temperature of the surrounding gas, provided 
that the opacity in the outer atmosphere is constant, 
i.e. K is independent of either gas pressure or temper- 
ature {a — (3 — 0). This is completely different from 
envelopes having outer convective region and occurs be- 
cause outer radiative zone decouples inner atmosphere 
from the surrounding nebula. As we demonstrated in 
>i5 5.11 this insensitivity of Mcr to the external nebular 
conditions holds also in a more general case of a = and 
arbitrary P (opacity independent on the gas density) , see 
equation (|53(l . This scaling is typical for the dust opacity 
which should dominate^ over the molecular opacity due 
to H2 and H2O under the conditions typical in the outer 
radiative zone, see opacity plots in Hayashi et al. (1979) 
and Mizuno (1980) and Figures El & H Thus, one would 
naturally expect envelope mass and critical core mass to 
be independent of po and Tq (except for the local value of 
opacity kq which may scale with local temperature Tq); 
this result is also completely insensitive to the structure 
of the innermost part (r < Ra) of the envelope, be it ra- 
diative or convective. Of course, as soon as a is nonzero, 
this degeneracy breaks and critical core mass starts to 
depend on Tq and po, as equation H54|) demonstrates. 

Previous studies self-consistently accounting for the 

^ Unless the dust opacity is very small, < 10^^ cm^ g"-*^ 



presence of the outer radiative zone (Hayashi et al. 1979; 
Mizuno 1980; Stevenson 1982; Nakazawa et al. 1985; 
Ikoma et al. 2000) have found Mcr ~ 10 Mq, depend- 
ing on the dust opacity and accretion luminosity of the 
core. These authors typically assumed constant M, i.e. 
Mc oc T, while we consider several possible accretion 
regimes taking into account the dependence of M on the 
core mass and its distance from the Sun. Note that for 
the same core formation timescale our assumed accretion 
law Mc oc (see Appendix yields higher M in the 
end of core accretion than constant M does; this acts to 
increase Alcr- If nucleated instability sets in when enve- 
lope mass is 30% of the core mass and kq = 0.1 cm^ g~^, 
as we assumed in all our numerical estimates, we find 
that critical core mass in the fast regime of planetesimal 
accretion is 53 Mq at 5 AU (present Jupiter's location) 
and 62 Mq at 10 AU (current Saturn's location). In the 
intermediate accretion regime Alcr is 30 M0 and 23 M0 
at 5 and 10 AU correspondingly; it drops to 1.8 M^ and 
0.7 M0 at these locations in the slow accretion regime. 
Ikoma et al. (2000) have found values of Mcr that are 
somewhat lower (sometimes by ~ 2) than those found 
in this work for the same M and kq because they used 
opacity in the outer atmosphere independent of P and 
T which has effect of increasing Menv (see below) while 
we use realistic dust opacity. Our neglect of atmospheric 
self-gravity is another reason for this difference. 

Inner envelope (r < Ra) can be either radiative or con- 
vective, depending on the detailed opacity behavior. We 
found that this part of atmosphere is radiative provided 
that a and /? in H13|l are such that the condition (|41|) is 
fulfilled; if Vq > Vad the inner envelope must be con- 
vective. Stevenson (1982) used constant opacity in his 
study which made the entire envelope convectively sta- 
ble and energy was transferred solely by radiation. On 
the other hand, in the important case of opacity dom- 
inated by small dust grains one would take a = and 

« 1 — 2. Equation l|41|) demonstrates that envelopes 
with such opacity should be convectively unstable if they 
consist of diatomic gas with 7 = 7/5, see Figures El & EI 

Of course, dust opacity cannot dominate in the whole 
envelope — at large depth pressure and temperature are 
so high that H2 and H2O opacities become more impor- 
tant, and their behavior cannot be described by a simple 
power law dependence ((T^ . Hayashi et al. (1979) and 
Mizuno (1980) showed that envelope is typically convec- 
tively unstable when these molecular opacities dominate. 
In addition to molecular k these authors also included 
constant dust opacity similar to Stevenson (1982), which 
had the effect of making atmosphere radiative in the 
outer region where k was dominated by dust and con- 
vective at greater depth where molecular opacity was 
more important. In their case this transition typically 
occurred quite deep in the envelope, way below what we 
would call Ra- Presence of such extensive radiative zone 
increases Menv and decreases Mcr- as long as a = /3 = 0, 
density in the radiative region varies as r~^ and atmo- 
spheric mass is evenly distributed in equal logarithmic 
intervals in r, see Appendix \n\ This augments total 
envelope mass compared to the mass contained in just 
the outermost part near Ra by additional factor equal 
to the logarithm of ratio of the outer to inner radii of 
this radiative zone (and this ratio is large). But if one 
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uses opacity with Vq > Vad then there is no radiative 
region below i?a, envelope becomes convective at much 
smaUer depth, right below Ra, and most of the atmo- 
spheric mass is concentrated near Ra, see Appendix IbI 
We expect the latter to occur [and Mcr to be higher than 
what Hayashi et al. (1979), Mizuno (1980) and Ikoma et 
al. (2000) have found] whenever opacity is dominated by 
small dust grains, since a = and /3 « 1 — 2 in this case. 
Convection at small depth in the case of dust opacity 
also wipes out the need to know the molecular opacities 
deep inside the atmosphere very accurately: behavior of 
K there is irrelevant as long as envelope is convective be- 
low Ra- 

6.2. Critical mass: implications for planet formation. 

Our results for the critical core mass necessary to trig- 
ger the nucleated instability fall within the range of 
previously estimated values of Mcr- In the terrestrial 
planet region, at 1 AU, we estimate^ Mcr ~ 20 M0 
if planetesimal accretion is in the intermediate or fast 
regime: RbX/R\ < 1 in both cases, meaning atmo- 
sphere with outer convective zone. Perri & Cameron 
(1974) and Wuchterl (1993), who studied convective en- 
velopes, would have obtained the same value of Mcr if 
they were to calculate it at 1 AU. 

At the same time, accretion in the intermediate or fast 
regime which keeps Rl high can only proceed for rather 
limited time span at ~ 1 AU until the isolation mass is 
reached (Lissauer 1993). Beyond this point protoplan- 
etary cores accrete basically at their geometric cross- 
section, i.e. in the slow accretion regime (Appendix 
Goldreich et al. 2004), for which Rb\/ R\ could be some- 
what higher than 1. After this transition occurs, critical 
core mass at 1 AU goes down to about 10 (see Figure 
Cj). Inward from 1 AU, in the region where "hot Jupiters" 
were discovered, Mcr decreases and reaches « 5 at 
0.1 AU. Since the amount of refractory material in the 
inner parts of protoplanetary disks is rather small and it 
takes very long time (~ 10* yr, see Chambers 2001) to 
collect this material into ~ 1 protoplanets, we con- 
clude that cores forming in situ in the terrestrial zone 
would not be able to undergo nucleated instability be- 
fore gas in the nebula has gone away (within ~ 10^ yr). 
These considerations strengthen the argument according 
to which "hot Jupiters" or their massive progenitor cores 
have migrated to their present locations from elsewhere. 

Although masses of existing terrestrial planets are 
clearly too low to drive the nucleated instability, they 
were high enough to retain quite substantial atmospheres 
while the nebular gas was still around. Our calculations 
demonstrate that Menv = 4 x lO^^ g, lO^'"^ g, 3 x 10^^ 
g, and 3 x 10^^ g for Mercury, Venus, Earth, and Mars 
correspondingly. This is to be compared with the present 
atmospheric masses of these planets: no atmosphere on 
Mercury, 5 x 10^^ g, 5.2 x 10^^ g, and 6.5 x 10^* g on 
Venus, Earth, and Mars. Clearly, the primaeval atmo- 
spheres of terrestrial planets have been heavily depleted. 
Massive primordial atmosphere can cause severe blanket- 
ing and melting of the core surface, effect which has been 

® Accurate value of Mcr in the terrestrial region can only be 
obtained after the effects of the vertical disk structure, differential 
shear, complicated flow geometry within the Ifill sphere, etc. are 
properly includ ed because in this part of protoplanetary disk Mc > 
Ml (see i|2 2.2l and Figure [7J. 



first considered by Hayashi et al. (1979). 

As we demonstrated in H6 6. II critical core mass in the 
region of giant planets sensitively depends on the accre- 
tion regime: at 10 AU it is less than M^ in the slow 
regime and as large as 60 M^ in the fast. At the same 
time, present mass of the Jupiter's solid core is estimated 
to be < 10 M^; current Saturn's core mass is between 
10 Me and 25 Mg, (Saumon & Guillot 2004). The ini- 
tial core masses of these two planets could have been 
higher — their total masses in high-Z elements can be 
as high as 30 — 40 M^ because some refractory materials 
can be dissolved in their envelopes. Transporting these 
elements from the core into the envelope would require 
some quite efficient mixing process such as vigorous con- 
vection (Stevenson 1982) and it is not clear at present if 
this is possible. 

Accumulation of the solid cores in the giant planet re- 
gion must have proceeded in the intermediate or fast 
accretion regime (or some combination of the two) be- 
cause smaller M would not allow solid cores to form be- 
fore the nebula dissipation this far from the Sun. This 
translates into values of Mcr which arc higher than the 
present core masses of Jupiter and Saturn. Thus, current 
data suggest that if cores were not initially more massive 
and subsequently eroded they were not large enough to 
trigger nucleated instability and accrete gas. The dis- 
crepancy is especially dramatic for Jupiter with its very 
low Mc- One way to resolve this paradox is to hypothe- 
size that nebular opacity kq was lower than 0.1 cm^ g~^ 
which brings Mcr down, see (I57II . Another possibility is 
suggested by the observation that core formation in the 
intermediate or fast regime can take rather short time, 
of order several Myrs, see (|A2(I . In this case, after the 
core of Jupiter has been mostly formed, there was still 
enough gas around to accumulate atmosphere. At this 
point, since the isolation mass has been reached, plan- 
etesimal accretion became considerably slower meaning 
that Mcr went down dramatically, below already accu- 
mulated Mc- As a result, gaseous envelope could not be 
supported by the energy release due to the significantly 
lowered M and it started to slowly contract on a thermal 
timescale. Thus, the transition between different plan- 
etesimal accretion regimes must be accompanied by the 
period of envelope adjustment similar to that described 
in Introduction. After Mcnv reached ~ Mc runaway gas 
accretion commenced. Similar scenario was studied nu- 
merically by Pollack et al. (1996) and Ikoma et al. (2000). 
We envisage the same picture to hold for Saturn as well, 
although its core must have taken longer to form than 
that of Jupiter and the nebula was appreciably depleted 
by that time; this may account for the smaller gaseous 
mass of Saturn. 

Assembly of the cores of ice giants — Uranus and Nep- 
tune — would have required even more time. Relatively 
small atmospheres of these planets containing only 1 — 4 
Me of H and He suggest that they either never under- 
went nucleated instability, or if they did this happened 
only after nebula has been very strongly depleted. In the 
former case all the gas that we see now in Uranus and 
Neptune has come from steady state atmospheres around 
the cores which were less massive than present cores of ice 
giants; after nebula went away these (sub-isolation mass) 
cores merged retaining their gaseous content (Genda & 
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Abe 2003, 2004) and producing Uranus and Neptune. In 
the latter case fast accretion leads to the accumulation of 
isolation mass (essentially the present day mass of ice gi- 
ants) and to subsequent nucleated instability before the 
nebula dispersal; however, if at that moment nebula was 
depleted by less than 10^ — 10'^, Neptune and Uranus 
would have much higher gaseous masses than they do 
now. 

Critical mass rather strongly depends on the envelope 
composition, namely on fi. Stevenson (1982) was the 
first to notice this fact for purely radiative atmospheres. 
Our equation H57(l confirms and generalizes this obser- 
vation — we find Mcr oc for any envelope with 
the outer radiative zone (independent of whether its in- 
terior is radiative or convective) . Dissolution of infalling 
planetesimals, erosion of the core by vigorous convection, 
or evaporation of some volatile icy content of the core 
can increase fj, and decrease Mcr considerably ("super- 
ganymedean puffballs", Stevenson 1984). On the other 
hand, envelope enrichment in high-Z elements might also 
increase opacity in the outer atmosphere and this can at 
least partially counteract the decrease of Mcr due to large 

6.3. Validity of assumptions. 

Because of the inherently analytical nature of this 
study aimed at singling out the most important aspects 
of atmospheric structure we have neglected a number of 
relevant phenomena which may be important in some 
cases. Among them are the dissolution of infalling plan- 
etesimals (Pollack et al. 1996), opacity jumps due to 
dust grain melting (Mizuno 1980), increase of planetesi- 
mal capture cross-section caused by the presence of the 
atmosphere (Inaba & Ikoma 2003), etc. We have also 
employed a set of simplifying assumptions such as hy- 
drostatic and thermal equilibrium of atmosphere, negli- 
gible gas accretion luminosity, and so on. In Appendix 
ICl we demonstrate these assumptions to be valid. Our 
treatment of convection relies on the absence of entropy 
gradient in the convective regions and in Appendix IdI we 
checked whether this assumption is appropriate. 

All our results are rather insensitive to the exact value 
of distance Rout at which atmosphere finally merges with 
the nebula as long as Rout ^ Rb- This is because atmo- 
spheric pressure outside i?s is not very different from Pq 
— planetary gravity cannot strongly perturb the pres- 
sure, while temperature is not very different from Tq: in 
the case of envelopes having outer convective zone devi- 
ation of T from Tq is directly related to the deviation of 
P from Pq by the condition of adiabaticity and, conse- 
quently, must be small outside of Bondi sphere. In the 
case of envelopes having outer radiative zone tempera- 
ture gradient is subadiabatic meaning that temperature 
deviation outside Rb is even smaller than in the convec- 
tive case. As a result, atmospheric pressure and tem- 
perature are not very different from their nebular values 
already at Rb and the exact location where P and T 
closely match Pq and Tq is not very important. 

Our use of opacity in the form H13(l may seem an 
oversimplification compared to other treatments (e.g. 
Mizuno 1980; Ikoma et al. 2000) which employ realis- 
tic opacity tables. However, as mentioned previously, we 
expect K in the outer part of the envelope to be domi- 
nated by dust and this typically (for /3 w 1 — 2) leads to 



convection inward from Ra obviating the need to know 
the gaseous opacity behavior at high temperatures and 
pressures. Some effect on the envelope mass and critical 
core mass might come from the change of equation of 
state caused by dissociation and ionization in the deep 
interior of the atmosphere, which e.g. might lead to the 
appearance of radiative regions at large depth. However, 
since atmospheric mass budget is dominated by the outer 
part of the envelope at r ~ Ra, change of k or 7 deeper 
down hopefully would not strongly affect the value of 

Menu • 

Thus, we expect our simple analytical treatment to 
provide robust qualitative picture of the envelope struc- 
ture and its dependence upon local conditions in the pro- 
toplanetary disk, and yield reasonable quantitative esti- 
mates of Menu and Mcr ■ 

7. CONCLUSIONS. 

We investigated steady-state structure of the atmo- 
sphere around protoplanetary core immersed in the 
gaseous disk. Our major results can be summarized as 
follows. 

Atmospheres split into those having outer convective 
or outer radiative (almost isothermal) zone. The former 
have entropy of the interior equal to the entropy of the 
surrounding nebular gas; the latter have interior entropy 
which is much lower than the nebular entropy, owing 
to the decoupling provided by the outer radiative region. 
Type of atmosphere around a given core is determined by 
the relationships between the Bondi radius Rb, photon 
mean free path A, and luminosity radius Rl- If envelope 
is optically thick at the Bondi radius {Rb ^ A) atmo- 
spheric type is set only by the value of the dimensionless 
parameter Rb\/R\, inverse of which has a meaning of 
radiative temperature gradient far from the core: outer 
envelope is convective when Rb\/R\ ^ 1, while when 
this parameter is > 1 atmosphere has an outer radiative 
zone. For the conditions typical in the protoplanetary 
disks (such as MMSN) cores having atmospheres of the 
first kind are common in the terrestrial planet region; in 
the region of giant planets atmospheres have outer radia- 
tive zone. Structure of the atmospheric interior is deter- 
mined by the dependence of opacity on gas temperature 
and pressure; it becomes especially simple (convective 
as soon as temperature starts to appreciably vary with 
depth) if opacity is dominated by small dust grains. 

In the terrestrial region critical core mass for nucle- 
ated instability Mcr depends only on the local values of 
nebular gas density and temperature. In the region of 
giant planets Mcr is insensitive to either po or when- 
ever opacity in the outer radiative zone is independent of 
the gas pressure, because outer radiative zone decouples 
inner parts of the envelope from the nebular gas (irre- 
spective of whether the atmospheric interior is radiative 
or convective); at the same time Mcr is a strong func- 
tion of accretion luminosity, opacity, and mean molecular 
weight. Critical core mass varies as a function of distance 
from the Sun because of the variation of po and Tq in 
the terrestrial planet region and because of the variation 
of planetesimal accretion luminosity and photon mean 
free path in the region of giant planets. Typical value 
of Mcr is several tens of M^ in the giant planet region 
(30 — 50 M0 at 5 AU) if planetesimal accretion was fast 
enough to account for the core formation prior to the gas 
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dissipation. Close to the Sun, at a ~ 0.1 AU, Mcr ~ 5 I am grateful for hospitality to Kavli Institute for The- 

independent of the planetesimal accretion rate. This oretical Physics where part of this work has been done, 
makes in situ formation of "hot Jupiters" by nucleated Useful discussions with Doug Lin, Peter Goldreich, and 
instability onto the locally formed protoplanetary cores Bruce Draine are thankfully acknowledged. Author is 
very unlikely. a Frank and Peggy Taplin Member at the IAS; he is 

Results of this study are also important for under- also supported by the W. M. Keck Foundation and NSF 
standing the ancient atmospheres of terrestrial planets, grants PHY-007092, PHY99-0794. 
which must have been significantly depleted. 

APPENDIX 

A. SUMMARY OF THE CORE ACCRETION RATES. 

Following Rafikov (2003b) we take the core accretion rate to be 

M = nUpRcRHO, (Al) 

where is a parameter set by a particular mode of accretion. In this work we consider the following three important 
accretion regimes. 

First one is characterized hy 9 k, p (see definition and assumes that core accretes planetesimals at a rate set 
by the core's geometric cross-section ~ R^ . This regime is valid when the random epicyclic velocities of planetesimals 
are larger than the escape speed from the core's surface and gravitational focusing is weak. We call this regime slow 
accretion. It may occur in planetesimal disks after cores have reached isolation mass by oligarchic growth (Chambers 
2001; Goldreich et al. 2004). 

Second regime of intermediate accretion takes place when random velocities of planetesimals are of the order of 
shear velocity across the Hill radius of the core Q.Rh- In this case gravitational focusing strongly increases accretion 
cross-section above its geometric value, and k, \. Note that this case corresponds to the boundary between the 
shear- and dispersion-dominated dynamical regimes (Stewart & Ida 2000); it also assumes vertical scaleheight of the 
planetesimal disk to be ~ Rh- This regime may occur during the oligarchic growth of protoplanetary embryos by 
accretion of large planetesimals (Kokubo & Ida 1998; Rafikov 2003b). 

Finally, third regime is realized when random velocities of planetesimals are so low that planetesimal disk is geomet- 
rically very thin and essentially two-dimensional. For that one needs random velocities to be smaller than p^^'^flRn, 
which leads to a rapid accretion with 9 « Such dynamically "cold" planetesimal populations can occur even 

in the presence of massive cores provided that some dissipative process such as gas drag can effectively damp plan- 
etesimal velocities. This can happen if fragmentation of large planetesimals in collisions at high velocities (induced by 
the gravity of protoplanetary cores) is capable of channeling a significant amount of mass initially locked up in large 
planetesimals into small bodies; see Rafikov (2003a, b) for details of this scenario. 

It is likely that each of these accretion regimes can take place during some stage of protoplanetary growth, e.g. start- 
ing with intermediate, switching to rapid (or some mixture of rapid and intermediate, depending on the fragmentation 
timescale, see Rafikov [2003b]), and, possibly, ending with a slow accretion phase. Typical accretion timescalc Tacc in 
each regime is 

M M / Af \ f ^ °^10' 

/ lUr, \ I ^ -^q7 ^2^^ intermediate, (A2) 



fast. 



Note that in all three regimes listed here protoplanetary growth proceeds as Mc oc r , where t is the time. 



B. CALCULATION OF THE ENVELOPE MASS. 



As we mentioned in H5 5.31 because gas density drops exponentially outside i?a, most of the atmospheric mass is 
confined within this radius and we may replace Rb with Ra in definition H51|l . Using H28|l and (|29|l for envelopes with 
radiative interior, and H42(l . I|43(l . and (|45(l for envelopes with convective interior we arrive at equation H52|l. in which 
^l/i is defined by 



^Sj,(u,w,O^Cw'^ (^^^' ^ j z^l^-^-l^ dz, (Bl) 





and for envelopes having radiative interior (Vq < Va^) 

/4V \ ^/(^+") 1 
C=f^j , u = t w = Vo, C = ;^-l, (B2) 

while in the case of envelopes having convective interior (Vq > Vad) 

/ vy \ 1/(4-/3)/. \l/(l+a) / \ 1/(4-/3) -, 
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Constant C is a power law ind ex of density dependence on 1/r, see equations (|28|) . (|29() and H44|l . 

Whenever C < 3 integral in (|B1|I is dominated by the contribution from r ^ Ra, i n wh ich case Me„„ depends only 
weakly on the lower integration limit Rc ^ Ra (this is why we set Rc/Rb to in |B1| ). Radiative envelopes with 
constant opacity (a = /3 = 0) having Vq = 1/4 and C = 3 contain equal amount of mass per every decade in radius; 
in this case formula (|52(l still describes the behavior of M^nv but with '^i ^ \n{Ra/Rc) replacing I|B1(I . This can 
considerably increase the envelope mass since Ra S> Rc ■ Condition C < 3 is always satisfied for dynamically stable 
convective envelopes which ought to have 7 > 4/3. In the radiative case one needs Vq > 1/4 for C < 3; whenever 
this is not fulfilled the envelope mass is dominated by the innermost part of the atmosphere near Rc, and Mcnv does 
depend on Rc- We do not consider radiative atmospheres having C > 3 deep in the envelope in this study. 



C. THERMAL TIMESCALE OF THE ENVELOPE. 

We calculate the thermal (or Kelvin-Helmholtz) time for the atmosphere with the outer radiative zone (see ® and 
convective interior (a = 0, /? = 1, Vq = 1/3 > Vad) with 7 — 7/5. We define thermal time as Tth = \Etot\/L, where 
Etot = Eth + Egr is the total energy contained within Ra — sum of the thermal and gravitational energy of the 
envelope. It is easy to verify that Eth ^ \Egr\, meaning that Etot ^ Egr as well. We calculate Egr using H42(l . (|43|l . 
and (|^ : 

/M 
G— X AT:p{r)r'^dr 

Ra 



;-^'^"«Hiej Uv^^tTJ t^A'-^) <'^'* 

This expression is valid for envelopes with convective interiors whenever 7 < 3/2, in which case gravitational energy 
bu dget is dominated by the innermost part of the envelope, near the core surface (see Harris 1978); this is why Egr 
in IjCip explicitly depends on Rc. This would be different for atmospheres with convective interiors having 7 > 3/2 
— th en the energy content is dominated by the outer parts of the envelope'' and is much smaller than that given by 
(fCTll . Ener gy is also small for envelopes having outer convective zone — the ir energy is low because of the high interior 
entropy and associated low density. In a sense, the specific estimate (|C1|) of Egr sets an upper limit on Tth and the 

degree of envelope's deviations from the steady state. 

Using our adopted values of 7, Vq, ^ad and luminosity @, HA1|I we find that 

5/3 /n 1 r.rr.2„-l\ \ ^ ^fo^'^, SlOW, 



Mc \ ' (0.1 cm^g 



Mm J \ Ko 



Tth ^ [ YT~ ) ( ~ — ) { yr aH^'^, intermediate, (C2) 



1 yr ajp^^, fast. 



Comparing this with (|A2|) one can see that typically Tth ^ ^acc justifying our quasi-static approximation to the 
treatment of the envelope structure. Massive embryos (Mc > 10 M^) accreting in the slow regime at a > 10 AU 
should have this condition violated and this may pertain to the development of nucleated instability in the region of 
giant planets, see H6 6.21 

At a given location in the nebula the energy stored in the envelope depends only on Mc. Since Mc changes due to 
planetesimal accretion, envelope energy should also vary in time giving rise to additional luminosity Lg caused by gas 
accretion: 

OT dMc Tacc 

Using the definition of Tth we then find that Lg/L ^ Tth /Tacc. Thus, whenever the quasi-stationary approximation 
(i.e. Tth ^ Tacc) is valid, gas accretion luminosity Lg is small compared to the core accretion luminosity L, and we can 
safely neglect it. Mass accretion rate of gas Mcnv ~ M {Mcnv / Mc) is always smaller than planetesimal mass accretion 
rate M if core is subcritical, i.e. M„„„ < Mr. 



D. EFFICIENCY OF CONVECTIVE TRANSPORT. 

Our use of equation relies on the assumption of convection so efficient that even infinitesimal deviation of 
temperature gradient from Vad is enough to transport the energy fiux produced at the core surface. If this is not the 
case one has to use mixing-length theory (Kippenhahn & Weigert 1990) to determine the value of V; here we check if 
this is ever necessary. Following Kippenhahn & Weigert (1990) we introduce 

Y7 Y7 W Y7 Y7 TT ad CjT^ 

X=V ~ Wad, W = Vrad - ^ ad, U ^ 



rf Pnp \GMcH^ 



1/2 



IQnaGMcT'^' dlnP' 



The latter is also true for envelopes with radiative interiors having Vo > 1/3 (which can exist only for 7 > 3/2); most of the energy in 
the radiative envelopes with Vq < 1/3 is near the core. 
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where 77 ~ 1 is a mixing length parameter, Hp is the pressure scaleheight, and a; is a deviation of temperature gradient 
from Vad, which has to be much smaller than unity for to apply. Value of x has to be obtained from the following 
equation (Kippcnhahn & Weigert 1990): 

[y/x + m-uy = ^U{W - X). (D2) 

It is clear from this equation that x \ whenever C/ <C 1, and UW <^ 1, since under these assumptions x ^ (UW)"^^"^. 
Thus, we look for conditions under which [/ <C 1, UW <C 1. 

In light of the results of tl3 3.31 fc Bl we describe the profiles of temperature, density, and pressure in the convective 
part of envelope by 

r = T,<i>, p = /,,$V(7-i), p = n$7/(7-i), $^i + v,,^f:^-:^y (D3) 

where T^, pb, Pb are values of temperature, density and pressure at the boundary of the convective zone r^. In the case 
of an atm osphere with the outer radiative zone studied in ti3 3.31 one has ^ Ra,Tb = Tconv, Pb = Pconv, Pb — Pconv, 
and (|D3(I reduces to 144|l . In the case of atmosphere with the outer convective zone rb — Rout, Tb ^ Tq, pb — po, Pb — Po 
and yields gp. 

One expects deviations of x from zero to be most pronounced in the outer, low density part of the atmosphere which 
might not be capable of efficient mixing. This corresponds to « Ra in the case of outer radiative zone, and rb ~ Rb 
in the case of outer convective zone, but in both cases <I>(rf,) « 1. In the atmospheres of first type <C i?sA) we 
find for a = 0, /3 = 1, 7 = 7/5 using ^ that 



slow, 

intermediate, (D4) 
fast. 

Also Vrad ^ ^ ad at r « Ra, meaning that W{Ra) ^ 1 and UW ^ U{Ra)- Thus, in the region of giant planets (a > 5 
AU) convection is efficient in the envelopes of protoplanetary cores accreting planetesimals in the slow or intermediate 
regime. In the case of fast accretion deviations of V from V ad at the level of 0.1 — 1 might be expected in the outermost 

parts of convective zone, at r « Ra- 

In the case of envelopes having outer convective zone ^ Rb^) one finds for U{Rb) expression similar to ljD4|l 
but without R'^/RbX terms. Also, V^ad ^ R\/Rb^ » 1 at r = Rb- As a result, 



slow, 

intermediate, 

(D5) 

fast. 

These estimates demonstrate that in protoplanetary atmospheres around the cores at a ~ 0.1 AU, which are likely 
to have outer convective zones, deviations from purely isentropic convection may only be important for fast and 
intermediate accretion regimes. 

These are the occurrences when the temperature gradient can go superadiabatic at the outer edge of t he c onvection 
zone and one would get a more accurate estimate of V in the convective regions by actually solving l|D2|l . e.g. see 
Papaloizou & Terquem (1999). We do not expect this possible superadiabaticity in some localized parts of atmosphere 
to strongly affect our estimates of Afgnt) and M^r since deeper down in the envelope density increases, convection 
becomes more efficient at transporting the energy, and entropy gradient goes to zero. 



^^"^"^ ~ coKopIGAL \Rb\ Rl ) 

VK) liir^rj 3 X 10 4^a,o / , 

I 0.2 a^o' , 



U{Rb) 



0.03 



0.1 cm^g 



15/4 



W{Rb)U{Rb)~0.5 



L 



Po4Ri 



Mr 



2/3 10- 



-1/2 



0.03 
0.4 ai 



1/2 
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